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A d-wave superconducting ground state for a doped Mott insulator is obtained. It is distinguished 
from a Gutzwiller-projected BCS superconductor by an explicit separation of Cooper pairing and 
resonating valence bond (RVB) pairing. Such a state satisfies the precise sign structure of the t-J 
model, just like that a BCS state satisfies the Fermi-Dirac statistics. This new class of wavefunctions 
can be intrinsically characterized and effectively manipulated by electron fractionalization with 
neutral spinons and "backflow" spinons forming a two-component RVB structure. While the former 
spinon is bosonic, originated from the superexchange correlation, the latter spinon is found to 
be fermionic, accompanying the hopping of bosonic holons. The low-lying emergent gauge fields 
associated with such a specific fractionalization are of mutual Chern-Simons type. Corresponding to 
this superconducting ground state, three types of elementary excitations are identified. Among them 
a Bogoliubov nodal quasiparticle is conventional, while the other two are neutral excitations of non- 
BCS type that play crucial roles in higher energy/temperaure regimes. Their unique experimental 
implications for the cuprates are briefly discussed. 

PACS numbers: 74.20.-z, 74.20.Mn, 74.72.-h 

I. INTRODUCTION 

An important issue in the study of high-T c cuprates concerns how superconductivity can arise in a doped Mott 
insulator [T]. The d-wave pairing symmetry has been usually attributed to the reason that the electrons avoid a strong 
local Coulomb repulsion. But in a doped Mott insulator the on-site Coulomb repulsion is so strong that the Hilbert 
space of the electrons is also drastically altered. It is thus no longer sufficient just to focus on the relevant attractive 
interaction as in the framework of the BCS theory. Rather a fundamental change in the underlying electronic structure 
should be taken into account before one can meaningfully address the issue of high-T c superconductivity. 

The simplest straightforward method of incorporating superconductivity with the "Mott physics" is to construct 
a Gutzwiller-projected BCS superconductor. As first envisaged by Anderson]]] in 1987, this class of state is always 
insulating at half-filling as all the double-occupancy states get projected out and the original Cooper pairs in the 
BCS wavefunction become neutralized, known as the spin-singlet RVB pairs which are "glued" by the superexchange 
coupling (the state is to be referred to as the fermionic RVB state below). Superconductivity arises only away from 
half-filling when the RVB pairs start to move in the presence of, say, empty sites in the hole doped case, which become 
partially charged Cooper pairs. Here the Cooper and RVB pairings are no longer explicitly distinguishable. 

The Gutzwiller-projected d-wave BCS state has been studied^ [3] intensively as a class of variational wavefunctions 
for the doped Mott insulator. It is also the basis for developing the so-called slave-boson approach]!], which is an 
electron fractionalization description with the Gutzwiller projection replaced by emergent gauge fluctuations around 
the spin-charge-separated saddle-points. It predicted^ |5] the presence of a high-temperature pseudogap phase over 
a finite doping regime and a superconducting dome at lower temperatures, which are both qualitatively consistent 
with the later experimental measurements in the cuprates pj [8]. 

However, the important long-range antiferromagnetic (AFM) correlations are notably missing in the Gutzwiller- 
projected BCS state at half-filling. In fact, it is a spin liquid with low-lying fermionic excitations [TJ [5J [3], which is in 
sharp contrast to a long-range AFM order with bosonic spin-wave excitations governed by the two-dimensional (2D) 
Heisenberg model. Of course one may argue here that once in the presence of some finite concentration of doped 
holes, the long-range AFM order or correlations will disappear anyway, and hence the Gutzwiller-projected BCS/spin 
liquid state could become stabilized eventually as a competitive ground statepQH]. 

But it remains to be a real challenge to understand how the long-range AFM order/correlations of an antiferromag- 
net/Mott insulator can be effectively destroyed by the motion of the doped holes, and how the associated energy, albeit 
only a small fraction in the total superexchange energy [5], gets turned off in order to gain the kinetic energy of the 
doped holes. The issue at the heart of a dope Mott insulator is the competition between the kinetic and superexchange 
energies - if the AFM correlations are weakened, like in the Gutzwiller-projected BCS state, it would be much more 
favorable to the hopping of the doped holes; with the enhanced AFM correlations at low doping, on the other hand, 
the kinetic energy of the doped holes will get strongly suppressed. The novelty inherent from such incompatibility 
and competition between the hopping and the superexchange processes is thus expected [TU] to be responsible for the 
unconventional nature of the superconducting transition as well as a complex pseudogap phenomenon over a wide 
temperature regime above T c in a doped antiferromagnet/Mott insulator. 
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Therefore, to properly accommodate such novelty, which may provide a basic understanding of the rich and mar- 
velous pseudogap properties in the cupratesjjj |8], both the RVB and Cooper channels should remain generally distin- 
guished even in the superconducting regime, which implies the necessity for one to go beyond the simple Gutzwiller- 
projected BCS state approach to adequately address the interplay between the magnetism and superconductivity. 

Another important property that the wavefunction should obey is an altered statistical sign rule: although the 
electrons always obey the Fermi-Dirac statistics, which dictates that two electrons of the same spin cannot stay at the 
same lattice site, new statistics will emerge in a doped Mott insulator where the no double occupancy constraint 
further enforces that two electrons of opposite spins cannot occupy the same site. In fact, the t-J model at half-filling 
is totally "bosonized" in the restricted Hilbert space where the usual fermion signs are completely diminished: e.g., 
the ground state only possesses the trivial Marshall signs [9] [12] which can be easily gauged away. The nontrivial 
signs start to reemerge in the doped case, induced by the hopping of the doped holes, which is precisely described by 
the so-called phase string effect [T3HT5] in the t-J model. The corresponding sign structure is actually independent of 
temperature, dimensionality, and is of statistical nature [2J \TE\, which eventually recovers the full Fermi statistical 
signs only at high doping in the dilute electron limit. Physically the phase string effect also provides an accurate 
mathematical description of the so-called "unrenormalizable phase shift" first emphasized by Anderson early on|16j. 
The latter is a total phase shift added up from all the electrons in the ground state, in response to adding/removing 
an electron into/from the system, due to the strong on-site Coulomb repulsion. Consequently, the irreparable phase 
string effect/unrenormalizable phase shift will make the Cooper pairing, associated with the doped holes, intrinsically 
distinguished from the neutral spin RVB pairing caused by superexchange, again suggesting the necessity to go beyond 
the Gutzwiller-projected BCS state description. 

A superconducting ground state distinct from the Gutzwiller-projected BCS state has been previously constructed 
by the present author and coauthors 17 with incorporating the above-mentioned sign structure. Such a ground state 
can naturally reduce to an insulating AFM state at half-filling, which well accounts for the long-range AFM order as 
well as the short-range spin-spin correlations with a highly accurate variational superexchange energy, known as the 
bosonic RVB statej9][l7]. In contrast to the aforementioned fermionic RVB state, the doped holes are quite unfavorable 
to hop in the bosonic RVB (neutral spin) background which only involves the spin pairing between different sublattices. 
It was then shown[T7] that the doped holes will force a fundamental change in the RVB structure in order to gain the 
kinetic energy at finite doping, characterized by emergent "backfiow spinons" accompanying the hole hopping[T7]. It 
is these "backfiow spinons" that will be associated with the Cooper pairs instead of the original bosonic RVB pairing. 
Consequently both the bosonic RVB and Cooper channels remain explicitly separated in the superconducting state 
[cf. Eq. (84) in Ref. [17]]. Such a wavefunction description has demonstrated a rich complexity in the pseudogap 
regime as resulting from the competition between the RVB and Cooper channels. 

In this paper, we show that an important simplification in this approach can be made by realizing that the afore- 
mentioned "backfiow spinons" are actually fermionic upon a closer reexamination of the sign structure. In the previous 
formulation, they are described in the bosonic representation[17|, which causes unnecessary complications because 
the extra fermionic statistical signs are mixed with the intrinsic phase string effect (cf. Sec. II A2). As the result, 
we obtain a greatly simplified self-consistent description of both the ground state and excitations for the doped Mott 
insulator. 

The key results are summarized in Sec. II. The general form of the superconducting ground state is presented in 
Sec. II A, which is distinguished from the Gutzwiller-projected BCS superconductor by a novel separation of Cooper 
and RVB pairings. It precisely satisfies the altered statistical sign rule of the t-J model in the restricted Hilbert 
space, which is of mutual semion type instead of the Fermi-Dirac one. In Sec. II B, it is shown that such a new 
class of wavefunctions can be intrinsically characterized by electron fractionalization, where neutral bosonic spinons 
and "backfiow" fermionic spinons together constitute a two-component RVB structure. The low-lying emergent 
gauge fields associated with such a specific fractionalization are of mutual Chern-Simons type, whose origin can 
be directly connected to the precise sign structure of the t-J model. The fractionalization formalism also makes 
the manipulation of the ground state significantly simplified as the constituent subsystems are more conventional, 
governed by the effective Hamiltonians presented in Sec. II B2. Finally, corresponding to this superconducting 
ground state, three distinctive elementary excitations are briefly discussed in Sec. II C. Among them a Bogoliubov 
nodal quasiparticle is conventional, while the other two are non-BCS like neutral excitations that play dominant roles 
in higher energy /temperature regimes, controlling the superconducting phase transition and other exotic properties 
different from a conventional d-wave BCS superconductor. 

In Sec. Ill, a microscopic justification of the present approach is presented in detail. It is based on a full bosonization 
formulation known as the phase string representation[T4"] of the t-J model, in which the whole nontrivial sign structure 
is explicitly captured by a topological (mutual statistical) gauge structure. Then we show that such a formalism under 
the no double occupancy constraint leads to the introduction of the two-component spinons in order to adequately 
describe the microscopically distinctive superexchange and hopping processes. Such a new exact formulation provides 
a precise starting point that naturally results in an clcctron-fractionalized description of the superconducting ground 
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state at a finite doping, and a highly accurate bosonic RVB description of the AFM correlation in the zero doping 
limit. An effective theory of the elementary excitations is also obtained within the same framework. Finally, the 
conclusion and perspective are presented in Sec. IV. 

II. KEY RESULTS 

We present this section to summarize the key equations/results of the present work, which basically addresses the 
issue how the ground state of a Heisenberg antiferromagnet/Mott insulator can be turned into a superconducting 
ground state by doping. 

A. Ground state ansatz 

For comparison, let us start with the Gutzwiller-projected BCS state ansatz, proposed[T] for the t-J Hamiltonian 
on a 2D square lattice, given by 



|*rvb) = PgM-BCS) 



(1) 



where |g?-BCS) denotes an ordinary d-wave BCS state and Pq is a Gutzwiller projection operator enforcing the 
following no double occupancy constraint 



(2) 



Because of Pq, the Cooper pairing in |d-BCS) reduces to the neutralized RVB pairing[T] at half-filling, whereas at 
finite doping the Cooper and RVB pairings are not explicitly distinguished, as schematically illustrated in Fig. 1. 
By contrast, the superconducting ground state obtained in the present work may be formally written as 



| RVB) 



(3) 



in which |RVB) denotes a neutral spin background that always remains half- filled as a Mott insulator, whereas the 
Cooper pairs associated with the doped holes are created by annihilating Nh electrons from |RVB), which are in 
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FIG. 1: (Color online) Schematic illustration of the Gutzwiller projected BCS state given in Eq. Q. (a) Half-filling: singlet 
electron (Cooper) pairs reduce to the neutral RVB pairs (blue-colored bonds); (b) Hole doping: the RVB pairs become partially 
charged Cooper pairs as they can hop to the hole sites (e.g., as indicated by yellow arrows). 
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singlet, d-wave pairing with an amplitude gij. Apparently such a superconducting state automatically satisfies the 
no double occupancy constraint without invoking the Gutzwiller projection Pq as in Eq. In particular, it is 

distinguished from Eq. by an explicit separation of the Cooper pairing and RVB pairing at finite doping. 

Schematically the RVB pairing in |RVB) and the Cooper pairing in |^g) are illustrated in Figs. 2(a), and 2(b) as 
well as 2(c), respectively. One easily sees the distinction between the neutral RVB pairing in Fig. 2(a) and the Cooper 
pairing in Figs. 2(b) and 2(c). In the latter, a pair of holes are involved, in which each hole will be generally associated 
with a spin via an RVB amplitude [the blue-colored bond in Fig. 2(b)], since the hole is created by annihilating a 
spin whose RVB partner already pre-exists in |RVB). Note that the holes are mobile here and thus their spin partners 
are also automatically changing with the hopping. Two spin partners associated with a Cooper pair in Fig. 2(b) can 
become RVB-paired again to further gain the superexchange energy, which results in the configuration shown in Fig. 
2(c) and serves the driving force for the Cooper pairing. In this sense, the superexchange/RVB pairing provides the 
ultimate pairing "glue" for superconductivity. 

With the Cooper and neutral RVB channels being explicitly differentiated in Eq. ([3]), generally some nontrivial 
phase shift effect will emerge, as introduced by A^. Specifically, A/j is given by 

A fc ee £ (« • • •<) Mh,hr- ^^) e -< Ail+ ^ + - + V) (4) 

where nf = 1— ^] c ] a c ia > denotes the hole occupation number at site I, and (ph is a bosonic wavefunction symmetric 
with regard to the hole coordinates {lh} — h, h, ■ ■ •, ^N h , which is generally present to ensure gauge invariance of the 

phase shift fields. Here the phase shifts |f2; h j , associated with the holes, will directly act on the "ghost" spin liquid 

state | RVB) to monitor the background spin correlations. In other words, A^ represents an "entanglement" between 
these two channels whose physical implications and mathematical definition are to be shown below. 

1. Superconducting phase coherence 
Due to the presence of Ah, injecting a hole into the ground state |\&G(Af/i)) will generally induce a phase shift by 

Cia^cW) ~ e lfl '\V G (N h + 1)). (5) 

Thus the wavefunction overlap between the bare hole state and the true ground state of Nh + 1 holes crucially depends 
on e lQi . 
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FIG. 2: (Color online) Schematic illustration of the present ground state (|3| or (201, which is structurally different from the 
Gutzwiller projected BCS state shown in Fig. [l] (a) The bosonic RVB pairs in jRVB): each pair only involves spin partners at 
opposite sublattice sites; (b) Doped holes are created by annihilating [as indicated by the dashed arrows which also represent 
the backflow spinons in the fractionalization formulation ( 20 1] the spins at the hole sites, and their RVB partners in |RVB) 



automatically become associated with the doped holes; (c) These partner spins associated with the holes can also form RVB 
pairs, which facilitates the Cooper pairing of the doped holes. Note that the important sign structure of the wavefunction, i.e., 
Ah in Eq. jjjl, is not directly shown here. 
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Furthermore, by noting that the Cooper pairing amplitude already pre-exists via g^ in Eq. ([3]), the superconducting 
off-diagonal-long-range-order (ODLRO) is essentially determined by (see below) 

(cft^) oc (RVB| e< W ')|RVB). (6) 

Hence the superconducting phase coherence and the coherence of a Landau (or more precisely, Bogoliubov) quasipar- 
ticle will be simultaneously realized. In other words, a "normal state" obtained by disordering the phase shift factor 
e lf2i will be intrinsically a non Fermi liquid with a vanishing quasiparticle weight. 
In the following we provide a simple proof of Eq. Q by taking 

ifh = constant (7) 

in Ah without loss of generality. Note that Eq. ^ then reduces to 

|* G ) oc Z)^|RVB) (8) 

with 

ij 

and 

Aj = e-^'+^Jcit^. (10) 



As shown in Appendix A, one has y-Dy = O(Nh) according to Eq. (8 1. Then, so long as giy (i'y) is a short-ranged 
function of \i — j\, an ODLRO can be identified in \^q): 



g ij (V ij )=0(5) (11) 
with 5 as the doping concentration (S = Nh/N, where N denotes the total number of lattice sites). Generally speaking, 



Eq. (11) represents that the Cooper pairing amplitude is formed (with the pairing symmetry determined by gij). 
The true superconducting ODLRO, (c^Cji), is thus indeed determined by the phase coherence condition in Eq. (raj, 
where Cli sensitively depends on the spin correlation in |RVB). 



2. Sign structure 

The phase shift £li is quantitatively given by 

e" <ft * = e"i (•!-*?) , (12) 

in which 

" " ' (13) 



and 

$9=53^(0. (14) 

where 8i(l) = Imln (zj — z{) (z t is the complex coordinate of site i), and n\ a denotes the spin occupation number 
(with index a) at site I, which always satisfies the single occupancy constraint 



5>t = l (15) 
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acting on the insulating spin state |RVB). 

Then each spin in |RVB) will contribute to a ±7r vortex via 3>f/2 in Eq. (12) with itself sitting at the vortex core. 
Vice versa, each doped hole will be perceived by the spins in |RVB) as introducing a 7r vortex, also via 3>f/2, with the 
hole sitting at the core. It implies that a doped hole and a neutral spin satisfy a "mutual semion statistics" as the 
phase shift Hi amounts to giving rise to ± 7r when one kind of species continuously circles around the other one once. 
Note that the single-valueness of Eq. (12) will be ensured by combining with $°/2. Thus the total phase shift added 



up in A/j represents a nontrivial entanglement between the doped holes and background spins, which will decide a 
"mutual semion statistics" sign structure in I^g) that is fundamentally different from that of a BCS state satisfying 
the Fermi-Dirac statistics. 

One may examine such a sign structure by a thinking experiment in which a hole in |\&g) goes through a closed 
loop c. At each step of nearest-neighbor moving of the hole, a singular phase or n is generated via a phase shift 
Oj in Ah depending on j" or J, spin that the hole "exchanges" with. Note that A/, also produces other phase shift 
contributed by other spins not "exchanged" with the hole, but their effect disappears after counting the total Berry's 
phase acquired by the closed-path motion of the hole. In the end, one finds 

|*g)^(-1)^ (c) |*g) (16) 

in which (c) only counts the total number of the 4,-spins that the hole has "exchanged" with along the loop c on a 
square lattice. 

The same sign factor (— l) iv h( c ) has been previously shown to be the precise sign acquired by the hopping of a 
doped hole through a closed loop c in the t-J model, i.e., the phase string effect [T3U 15] , This effect is proven to 



be dynamically irreparable and is thus of statistics nature. So the phase shift f2 defined in Eq. (12 1 is necessarily 
generated by the motion of a doped hole, representing an emergent new statistics |llj in doped Mott insulators with 
the Hubert space restricted by no double occupancy constraint. 

It is noted that at finite doping, the exact topological sign structure identified based on the t-J model is generally 
given bv[T5] 

r c = (-i)W)t<«0 x (-l)^M (17) 

which appears, say, in the partition function 

Z = Y j t c Z{c) (18) 

c 

whe re Z(c) > for any closed path c of the multi-ho le / s pin configurations at arbitrary temperature. Compared to Eq. 
(16), an extra sign factor (— l) N h( c ) appears in Eq. (17) in which N^(c) counts the number of hole-hole exchanges on 

the path c. It is straightforward to verify that that fermionic signs of the doped holes created by ( 9ij c i^ c ji ) m 

Eq. ([3| can precisely account for such a sign factor. Consequently combined with the phase shift in A^, the nontrivial 
sign structure of the t-J model is naturally satisfied by the ground state (vPg) m Eq. (|3| so long as the neutral spin 
background |RVB) does not contribute to additional statistical signs as shown below. 



3. |RVB) as a spin liquid state 

As already mentioned, |RVB) describes a "ghost" spin state, which remains one spin at each site of a square lattice 
even in the doped case. But the spin state can evolve from an AFM long-range ordered one to a spin liquid with only 
short-ranged AFM correlations as the doping concentration is increased. 

Generically |RVB) can be expressed by 

|RVB) = ]T $rvb fa, oi, • • ;o N )c\ ai cl a2 ■ ■ ■ 4jO) (19) 
{*.} 

in the electron c-operator representation. In nature it is a bosonic state with the wavefunction ^rvb ({c s }) = 
^partition ll(ij) ("l) 1 ^^ f° r eacn given spin configuration {u s } = o\,o\,- ■ -,(7^. Here the RVB pairing amplitude 
(— lyWij connects two antiparallel spins denoted by i (up spin) and j (down spin), with the summation running over 
all possible pairing partitions for the given {a s }. The staggered sign (— l) 1 (the Marshall sign) is explicitly separated 
from Wij such that the latter remains a smooth function of the distance between even and odd lattice sites of a square 
lattice at different doping concentrations. 
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Such |RVB) is a generalized Liang-Docout-Anderson type bosonic RVB state, which can naturally recover the 
correct antiferromagnetism in the zero doping limit [9l fT7] . But a long-range RVB pairing in the antiferromagnetic 
phase will generally destroy the phase coherence condition in Eq. ^ because the ±7r vortices carried by the two 



spinon partners of a long-range RVB pair, according to Eq. (12 1, do not compensate each other and result in a 
phase disordering. Only can a short-ranged RVB pairing lead to a vortex-antivortex binding in Eq. ^ and thus the 
superconducting phase coherence. In other words, |RVB) has to become a spin liquid in the superconducting phase. 
How the RVB amplitude W^ evolves with doping and self-consistently becomes short-ranged in the superconducting 
state will be shown below. 

B. Electron fractionalization 

What is the physical implication of the explicit separation of the Cooper and RVB pairings in the ground state (J3|? 
In the following we show that it actually corresponds to a unique electron fractionalization. 

1. Ground state in electron fractionalization form 
The ground state |^g) in Eq. Q can be reformulated as a direct product state 

|* G ) =P(|$ h )®|$„)®|* 6 )) • (20) 

The coefficients, ifh, Qij, and Wij, appearing in the original |^g) (cf. Sec. II A), are incorporated into three subsystem 
states as follows: 



1> '2, 



■K h l 



■\0)h 



(21) 



and 



|$„) = exp ( }2& 3 al ia ] t 



|0)a , 



(22) 



as well as 



|$ fe )=exp[^W^ t 6^ | |0> 6 



(23) 



Here the bosonic wavefunction tp^ in |$^) defines a "holon" state with a bosonic creation operator h\ acting on a 
vacuum \0)h', \&b) defines a neutral "spinon" state with an RVB pairing amplitude Wij, where b ia as a bosonic creation 
operator acts on a vacuum |0)&; and |$ a ) defines a "backflow spinon" state with the pairing amplitude gij = (— l) 1 g%j, 
where a\ a denotes a fermionic creation operator acting on a vacuum |0) o . 

The projection operator P in Eq. ( 20 ) is defined by 



P = PrP* 



(24) 



in which P s will enforce the single-occupancy constraint Eq. (15) in the spinon state such that 

|RVB)=P S |$ 6 ) (25) 

with n b ia = b\ a bi a ; and Pb will further enforce 

n a lS = n^nl , (26) 



such that each a-spinon always coincides with a holon as J2<y n is 



according to Eqs. (26 1 and (151 (here 



- J 



is and 



,h - tf 



h]hi with a = — a). By applying P, the physical Hilbert space is restored in Eq. (20) as 



schematically shown in Fig. [2] in which the a-spinons are indicated by the dashed arrows at the hole sites, while there 
is always a 6-spinon indicated by a solid arrow at each lattice site. 

Note that the phase shift factor e~ lf2i in A/j has totally disappeared in the above direct-product expression Eq 
where the electrons break up into the fractionalized building blocks: the holon h\ the spinon b\ 

conventional 1 



(201, 



and the backflow 

The fractionalization form of the ground state Eq. (20) 



spinon a\ : forming three rather "conventional" sub-states 

can be straightforwardly obtained by substituting into Eq. Q the following decomposition form of the electron 
annihilation operator: 



= Ph. 



with 



htaU-afe^ 



(27) 



(28) 



which acts on the insulating "vacuum" \0)h <£> |0) a (8 |RVB). On the other hand, in the neutral spin state |RVB), the 



-operator can be reexpressed in terms of the bosonic spinon operator b ia to result [17] in Eq. (25) from Eq. (19), 



according to the decomposition Eq. (48) given in Sec. Ill, where it is further demonstrated that thefull spin operator 
can be expressed as 



PS, 



(29) 



with 



Si — Sj + Sj , 



(30) 



where denotes the spin operators for &-spinons [defined in Eqs. (52) and (53)] and Sf for a-spinons [defined in Eqs. 
((73} and iffl'. 



2. Effective Hamiltonian 



The electron fractionalization form [Eq. (20)] of the ground state I^g) wm make the manipulation of the super- 
conducting state more easily than in the original form [Eq. ([3])]. 
Define 



|* G ) = |$0 <g> |$b 



(31) 



such that |* G ) = P\^g)- Then c icr |*G) = -Pc i(T |* G ) and Sj|* G ) = -PSi|* G ), in which c ia and Si directly act on the 
fractionalized states. 

Based on the t-J model, we find that the direct product state I^g) in Eq- (31 ) can be effectively determined as the 
ground state of the following effective Hamiltonian 



H c ff — Hh + H s + H a 



(32) 



which is composed of a holon hopping term 



Hh — —t. 



(33) 



a 6-spinon pairing term 



(34) 
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and an a-spinon term 



H„ — —t„ 



{ij) 



+jJ2&-s b J + s b l 

{ij) 



(35) 



where Ay = J2 a aa ia a j~<j ana - Vij = ~K — ) f° r 3 = is a d-wave sign factor. Note that the chemical (Lagrangian 

multiplier) terms implementing 



E 



(36) 



and 



(37) 



are all omitted in Eqs. (33), (34), and (35) for simplicity. One can always add them back in real calculations. 
Based on Eq. (31 1, the parameters, th ~ t, J s ~ J, t n ~ and J Q ~ J 



A" 



in Eqs. ( 33 )-( 35 1 can be determined 



as variational parameters minimizing the ground state energy of P\^q) at a given doping concentration, which will 
involve the projection P and whose detailed magnitudes will not affect the general consequences to be outlined in the 
next subsection. 

In Hh and H s , the /i-holons and 6-spinons are generally coupled to the U(1)®U(1) gauge fields, A|- and Ah, 
respectively, in Eqs. (33) and (34), which are topological (mutual Chern-Simons) fields as their gauge-invariant flux 



strengths in an arbitrary counter-clockwise closed loop c are constrained to the numbers of spinon and holon matter 
fields within the enclosed area E c , respectively, 



E^ 



'it 



and 



E4 



/es c 



(38) 



(39) 



The link variables, Afj and Ay, can be regarded as mediating the mutual statistics coupling between the charge and 
spin degrees of freedom, i.e., the "mutual semion statistics" entanglement introduced by the phase shift factor e~ ll! 



in the original ground state Eq. (|3|). In addition, the constant link field 
flux per plaquette, which is originated from $^ term in Eq. (12). 



in Eq. ( 35 ) describes a non-dynamic tt 



3. Ground state as a mean-field solution 



Then the ground state (31 ) as a self-consistent mean-field solution of H e & in Eq. (32) can be constructed as follows. 
First of all, suppose the holon state \<&h) governed by Hh in Eq. (33) become Bose-condensed [cf. Eq. Q]. Such 



holon condensation will then lead to Ah 



terms of Eq. (39) 



Ah, with Ay depicting a uniform flux ^ D Ah = irS per plaquette, in 



= if both i and j be 



Then H s in Eq. (34) can be diagonalized, resulting in a mean-field solution given in Eq. (23), in which W t 



ong to the same sublattice and decays exponentially at large spatial separations for opposite 
r |17j. Here ry is the spatial distance and £ is the characteristic pair size 



sublattice sites i and j: |Wy| oc e 

determined by the doping concentration: £ = ay ^ (a is the lattice constant). Hence, the spin background |$&) 

indeed becomes short-ranged at finite doping with a finite spin gap E g oc (5J|10| 117] . Once the 6-spinons are all 
short-range paired up in the fluctuations of A|- would become negligible for the long- wavelength physics, i.e., 
Y^c^tj ~ ^ ^ 0r a large loop c as compared to £, according to Eq. (38). Self-consistently, the two subsystems of the 
holons and 6-spinons are decoupled as depicted by |$&), as the ground state of 



Hi, + H s 



(40) 
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which is known as the phase string model[10, 18 or mutual Chern-Simons gauge theory model 19 . 

It is interesting to point out that |RVB) = P s \$b) is of the same form as the Liang-Docout- Anderson type RVB 
wavefunction at half-filling, which has been previously proposed [9] as a very accurate variational ground-state wave- 
function for the Heisenberg model. Indeed, at half-filling, in the absence of holes, |^g) simply reduces to |RVB), 
with £ — > oo or obeying the power law at large spatial separation of ij: |Wj,-| oc 1/ | r^- 1 3 [T7| . The ground state 
energy and staggered magnetization of the antiferromagnetic ordering determined numerically based on such |RVB), 
obtained from H s , are highly accurate as compared to the exact numerical results [TJ]. This indicates that the bosonic 
RVB mean-field description in Eq. (34), which reduces to the Schwinger-boson mean-field theory at half-filling 20J , 
has accurately captured both short-range and long- wavelength correlations of the Heisenberg model in this limit. 

Finally, note that the fermionic backflow a-spinons in H a [Eq. (35)] are gauge neutral without coupling to the 
internal mutual Chern-Simons gauge fields. Here the last scattering term with the 6-spinons in Eq. ( 35 1 can be safely 
omitted in determining the ground state |$ a ), due to the above-mentioned gap E g opening up in the spin excitation 
involving 6-spinons in Then the bili near terms in Eq. (35) can be straightforwardly diagonalizcd with a proper 

gauge choice of </>^-, to result in Eq. (22) with a d-wave amplitude g^j. Due to the presence of contributing to 
a 7r-flux per plaquette, the a-spinons will form Fermi pockets at both (0,0) and antinodal point (7r,0), etc., and a 
staggered current loop is expected to be present at A°- ^ 0. The physical implications of the a-spinon excitations 
will be further discussed later. 

To end this section, let us examine the electron pairing parameter Afp = (cifCji) based on Eq. 
expressed in the present fractionalized state by 



(28). It can be 



A 



so 



(41) 



Namely, the pairing amplitude and symmetry will be determined by the pairing order parameter of the a-spinons, and 
the superconducting phase coherence is decided by Eq. (|6]). The latter is realized as the RVB pairing of the 6-spinons 
in |$b) becomes short-ranged with a finite £, such that the n- vortices and -antivortices attached to them, according 
to Eq. (13), are all confined to form the vortex-antivortex pairs. Namely, superconductivity will be protected by 
a "ghost" spin liquid state. As noted before, the superconducting phase coherence can either disappear as £ — > oo 
either in the long-range AFM state near half-filling, or in the overdoped regime when the RVB pairing in |RVB) is 
diminished by doping. 



C. Elementary excitations 



Once the ground state ansatz Eq. ^ or its fractionalization form Eq. ( 20 ) is determined, the corresponding 
low-lying elementary excitations, which r eflec t the novel correlations in the ground state, will also nat urally manifest. 

In fact, the effective Hamiltonian Eq. (32) determines not only the ground state |^g) m Eq. (31), but also some 
nontrivial excited states. In the following we first show the existence of two novel elementary excitations which are 
uniquely governed by H e g. Then we show that a conventional Bogoliubov quasiparticle excitation will also appear as 
a collective mode that goes beyond H c g , which will remain protected within the basic characteristic energy scale E g 
in this non-BCS superconducting state. 



Spin-roton excitations 



In the ground state Eq. ( 20 1 , one does not see the trace of the electrons directly - such a strongly correlated electron 
system seems entirely fractionalized, as described by the bosonic RVB paired spinons, Bose condensed holons, and 
d-wave paired backflow spinons, which form a direct product (a generalized "spin-charge separation") state. 

However, we find that the single 6-spinons and holons will not be truly present in the low-lying energy spectrum 
to become real elementary excitations. This is because 6-spinons and holons are not gauge neutral - they carry the 
"gauge charges" of the mutual Chern-Simons fields, A h and A s , while provide the "topological sources" to generate 
A s and A h , respectively, and thus their excitations, by breaking up the correlated patterns formed in the ground state, 
will generally invite nonlocal responses from the whole system, which would make such excitations too costly [2"T]. 

For instance, one can imagine spinon excitations created by breaking up an RVB pair in |$b), described by the 



effective Hamiltonian H s in Eq. (34). However, each unpaired 6-spinon will induce vortexlike superfluid currents via 



A s from the condensed holons according to Hh in Eq. (33), leading to the so-called spinon-vortex composite object 



which is logarithmically divergent in energy, as discussed in Refs. |22H24j . Therefore, these spinons can only exist 
in the RVB pair condensate in the ground state, where such vortex currents get effectively cancelled out due to Eq. 
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( [38] ) , but not as a single excitation at low energy. Namely, single o-spinon excitations will be "confined" [TO] in the 
bulk of the superconductor. 

On the other hand, an integer (S = and 1) spin excitation involving a bound pair of spinons excited in are 
still allowed [Ml 125] . in which the effect of vortex and antivortex bound to individual spinons get cancelled out in the 
long distance in Hh such that its excitation energy becomes finite with a mean- field gap E g according to H s . 

Such a neutral spin mode carrying an either integer S — or 1 quantum number, is called a spin-roton|25j. The 
spin-roton excitations will not destroy the phase coherence at finite temperature until T c , where the spin-rotons 
disassociate into free spinon-vortices '22-24]. It has been shown [2~5] that a simple T c formula: 

T c ^§- (42) 

can be determined with E g ~ 8J denoting the core energy of the spin-rotons, degenerate for S = and 1, which 
is in excellent agreement with the experiments, with the "resonancelike" modes observed in the Raman A\ g chan- 
nel and neutron scattering measurements consistently interpreted as the spin-roton excitations with S = and 
1, respectively 26 . It also provides a natural explanation why the two modes are energetically degenerate in the 
experiment |26) . 

Here the spin-rotons are the most essential elementary excitations of non-BCS-type above the superconducting 
ground state Eq. ([3] or Eq. (20), which directly controls the superconducting phase coherence Eq. ^ via a 



characteristic energy E g , In the AFM long-range ordered state near half-filling, one has E g — > such that T c 
vanishes, and the S = 1 spin-roton excitation will naturally reduce to the gapless spin wave. 



2. Fermionic a-spinon excitation 



In the ground state Eq. (20), there are two distinct branches of spinons. The 5=1 spin-roton excitations related 
to the dynamic correlation function of have been discussed above. According to Eq. (30), the backflow a-spinons 



will contribute to another branch of S = 1 excitations as governed by H a in Eq. ( 35 ) . Here the single a-spinons are 
gauge- neutral and can be excited by breaking up the d-wave pairs in |$ a ) [Eq. (22)], which is expected to provide 



a characteristically different spectral contribution below the spin-roton "resonance" modes mentioned above [at E > 
Eg, such an S = 1 mode composed of the a-spinons may strongly decay into the spin-roton mode via the last term in 
Eq. p5l]. 



Besides contributing to the spin spectral function, a single a-spinon can also directly appear in the single-particle 
channel. According to the decomposition Eq. (28), a coherent term may emerge as the first term in 



-Ka\ s {~aye^+:h\:a\ s {-a) 



(43) 



with the holon condensation (hi) = hg (: h\ := h\ — Hq) and the superconducting phase coherence (e lSli ) ^ 0. In 



other words, the a-spinon excitation may be directly probed by ARPES as a coherent term appearing below T c with a 
weight \h \ 2 cx 8 and disappearing above T c when ^e l ^ ; ^ = 0. Such a coherent term will be nonetheless distinguished 

from the true quasiparticle excitation, which will be obtained as a coll ecti ve mode, i.e., a bound state of the holon and 
a-spinon by the singular phase shift fli from the second term in Eq 
Since the low-lying a-spinon excitation will appear near the antinoda 
the true quasiparticle, after absorbing the phase shift field e* n 



43), as to be given in the following subsection. 



region [i.e., the momentum (0,7r), etc.], while 
is a nodal quasiparticle around (7r/2, tt/2), there is 
a "dichotomy" between these two kinds of excitations in the single-particle channel, which will be explored in detail 
elsewhere. 

Finally, it is noted that in a strong magnetic field, the d-wave pairing of the a-spinons may be first broken down 
by the Zeeman energy before the spin-rotons (of energy > E g ) in |RVB) get excited to destroy the phase coherence 
Eq. ([61. In this case, the a-spinons in Eq. |35l will form coherent Fermi pockets with A° = such that the 



superconaucting order parameter in Eq. (41) can also vanish. Then a new normal state characterized by small Fermi 



pockets of the a-spinons can be realized by applying a sufficiently strong magnetic field without encountering the 
phase disordering boundary. Of course, such a T = transition caused by the Zeeman effect should be compared 
to another possible route to a normal state via the generation of new type of magnetic vortices at strong magnetic 
fields!' 
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3. Quasiparticle as a collective mode 
So far we have discussed two novel elementary excitations in the superconducting state, i.e., the spin-rotons and the 



a-spinons, which are apparently non-BCS-likc as determined by iJ ff in Eq. (32). In the following we point out that 
the conventional Bogoliubov quasiparticle excitation will reemerge as a collective mode in the full t-J Hamiltonian, 
although it is not an eigen solution of H c g. In other words, a Bogoliubov quasiparticle can be regarded as a bound 
state of the fractionalized building blocks, which are glued by the residual interaction in the original t-J model [2T | 127 ) . 
However, such a quasiparticle excitation will be stable only in a sufficiently long-range, low-energy regime, where the 
superconducting state will still behave like a conventional d-wave BCS superconductor as other exotic modes are not 
yet excited. 

In order to generally trace a quasiparticle excitation, one may directly create a bare hole (particle) by the electron 
c-operator on the ground state I^g)) an d then follow its behavior via the following equation-of-motion (cf. Sec. Ill 
D for details): 

— i9 t Ckcr — — (ek — / i ) c k<r ~~ Ak<rcL k + decay term + scattering term (44) 
which is obtained after a linearization in terms of the mean-field order parameters. Here Ak cx A k is given by Eq. 



(94). The decay term is given in Eq. (95) which corresponds to the process discussed in the above subsection that a 



doped hole dissolves into an a-spinon, i.e., 



decay term ~ ( e lfis ) /iS a Lr (45) 



in the superconducting background of Iiq =^ and (e tVti j =^ 0. Such a coherent mode will appear in the antinodal 
region where the Fermi pockets of the a-spinons locate, and disappear once the superconducting phase coherence gets 



0. The scattering term reads [cf. Eq. (96) 



scattering term = ^ (2ir k+q - Jr q ) [ac k+cla S^ + c k+q _ ff S q - ff ] (46) 
q 

with Tk = cos k x a + cos k y a, which represents a process that the doped hole scatters with the spin-rotons excitations 
created by above the resonancelike energy E g . 

Therefore, if we only focus on the low-lying excitation below the spin-roton energy E g and around the nodal region, 



the scattering and decay terms in Eq. (44) can be all neglected. Then Eq. (44) and its quasiparticle counterpart can 
be combined tc 
which leads to 



be combined to give rise to an elementary excitation of the Bogoliubov quasiparticle type a ko . cx u\ L c\ i(J + crvkC-k-a, 



id t atj* G ) = E k atj* G ) (47) 



and ako-l^c) = 0, with itk, «k> and the energy spectrum E^ = y (e^ — fi) 2 + (Ak) given in Sec. Ill D, not different 
from a usual d-wave nodal quasiparticle in a BCS framework. Here the chemical potential fi is determined by requiring 
Skc c ko- Ckff = (1 — and one arrives at a very important conclusion that a d-wave Bogoliubov quasiparticle 

excitation is still well preserved in the present non-BCS-like superconducting state. It is not given by the effective 
Hamiltonian H Q g, but emerges as a collective mode of the original t-J model, which is ensured by the superconducting 



ODLRO (41) protected by a finite minimal spin-roton energy E 



III. MICROSCOPIC JUSTIFICATION 



In this section, we justify the ground state ansatz, Eqs. (|3| and (20), as well as the elementary excitations outlined 
in Sec. II, based on the t-J model. We shall start with an exact reformulation of the t-J model in an all-boson 
formalism in which the hidden sign structure can be explicitly revealed. Such a precise sign structure will then play 
an essential role in determining the peculiar structure of the ground state and elementary excitations in terms of the 
electron fractionalization. 
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A. Phase string representation of the t-J model 



The t-J Hamiltonian is a minimal model of doped Mott insulators, with the Hilbert space restricted by the no 
double occupancy constraint ^ in the hole-doped case. It has been rigorously demonstrated that the fermion signs 
are completely diminished at half-filling due to the "Mottness" enforced by the constraint where the t-J model 
reduces to the AFM Heisenberg model. The residual fermion signs only start to re- emerge upon doping, which can 
be mathematically described by the so-called phase string effect [T3TH5] . 

In order to explicitly keep track of such a sign structure, the phase string representation of the t-J model has 
been previously introduced p3], in which the electron annihilation operator can be fully bosonized by the following 
decomposition 

da = h\b la e 1 ^ (48) 

in terms of the bosonic holon creation operator h\ and the bosonic spinon annihilation operator b la , with the phase 
string effect explicitly embedded in the phase factor 

e**i~ = ( _ ff )V* [*«■-*?-"*?] (49) 



In e lQia , $f and $° are defined in Eqs. (13 1 and (14), respectively, and is given as follows 

*« fc =$>aK (so) 

which is nonlocally associated with the holon occupation number at site I. Here e lQia also plays a role to restore 
the fermionic statistics of c ir7 , in the Hilbert space restricted by the single occupancy constraint 

5>l +«? = !• (51) 

(7 

Correspondingly the spin operators = S b t with are defined by 

S(* = ^X> 6 U" ( 52 ) 

and 

S\+ = (-l)%b iie ^, St = (-l)% bif e-^ (53) 

which involve the holon degree of freedom via $^ defined above. 
Then the t-J Hamiltonian 



flt-J = H t + H, 



J 



= -t E 4 a c ja + h.c. +jJ2(s % -S 

(ij)<r (ij) 

under the constraint n, = J~) c\ a Ci a < 1 can be reformulated as|14j 



H t = "* E { h l h ^ Ah ) (blh^* +laA >') + h.c. (54) 



and 



J 2 

(ij) 



E^i)^' ( 55 ) 



with the bosonic RVB order operator 



x E ( " ? r h <^ » ( 56 ) 
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for the NN sites. 

Here the unique feature is the presence of three link fields: Ay, Ay, and (py, which capture all the nontrivial 
signs inherent to the model (without them, the model would simply reduce to a pure bosonic one without any phase 
frustration). They are defined by 



A1 



E 



m) 



an 



la 



(57) 



and 



(58) 



(59) 



transformations 



The flux strengths of Ay and Ay are given in Eqs. (38) and (39 1, respectively, which are invariant under the gauge 



and 



hie 1 



birr -> hi, 



A" 



-)• A'y + fa - tpj) 



4 h 



Ah 



(60) 



(61) 



Thus the t-J Hamiltonian in the phase string representation has an intrinsic U(1)<S>U(1) gauge structure and according 
to Eqs. (60) and (61 1, the holons and spinons will carry the gauge charges of the gauge fields, Ay and Ay, respectively, 
the strengths of Af 



and 



On the other hand 
respectively according to Eqs 
lattice without any dynamics. 

Therefore, the full sign structure hidden in th e t-J mod el h as been explicitly sorted out and captured by the link 



are generated by the local densities of the spinons and holons, 
( 38 ) and ( 39 ) . Finally simply describes a uniform 7r flux per plaquette on a square 



variables, Ay and Ay, as well as <j>y, in Eqs. (54) and (55) of the phase string/bosonization representation 
half-filling, Ay vanishes such that there is no nontrivial sign left in Hj, while H t = under the constraint (51) 
the fermionic signs of the electrons totally disappear here, and the residual intrinsic signs only re-emerge as t 



At 
So 

,he holes 

are doped into the system, which are precisely represented by the aforementioned topological gauge fields, and in this 
sense the t-J model becomes an intrinsic gauge model in the phase string representation. 



B. Electron fractionalization 



Hence the phase string representation of the t-J model constitutes a suitable starting point, as it smoothly connects 
the doping problem with the undoped antiferromagnet - the latter can be well described by the bosonic RVB state 
without involving any sign problem [51 117], 

However, as pointed out in Ref. |17) . a bosonic spinon defined in the phase string representation, satisfying the 
constraint (51), is not strictly charge-neutral: i.e., it is involved not only in the "neutral" superexchange process 
described by Hj [Eq. (55)], but also in a "backfiow" process accompanying holon (charge) hopping in the H t term 
[Eq. (54)]. As a matter of fact, it is shown [T7] that the bosonic RVB pairing is incompatible with the backfiow 



accompanying the holon hopping (see Sec. Ill C below), and thus it is important to further distinguish these two 
different types of processes for spinons in this particular representation. 



Two-component spinon description 



In order to properly accommodate these two distinct correlations of spins, a neu tral spinon Hilbert space has been 
introduced|17) previously, which satisfies the constraint n\ a — 1 instead of Eq. (51 ), even in the doped case. Such 
a "ghost" neutral spin state, denoted by |RVB), can form a direct product state with the holon state |$^) 



|RVB) <g> |$ h ). 



(62) 
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Then the excessive spinons coinciding with the holon sites in Eq. ( 62 1 should be removed as they violate the Mott 
constraint (51 ). A physical state constrained by Eq. (I5T|) may be then realized by 



|*) = n B |RVB)® |$ h ) 



(63) 



annihilating the unphysical neutral spinons at the holon sites. Here I1b 



With ft B = • • ■ {bl s a s ni s )( b h + l<?s + i n l s + 1 ) 

may be regarded as introducing a new type of spinons, namely, the backflow spinons [T7]. in the "vacuum" state (62 1. 
Below we show how to mathematically accurately incorporate such two-component spinon building blocks into the 
phase string representation. 

Define as the creation operator for a backflow spinon, via a one-to-one mapping 



(fy CT nf) |RVB) ® |$ h ) m. P B 



(aLe^) |0) o ® |RVB) ® |$ h ) 



(64) 



where |0) a denotes the vacuum state of the et-spinons and the projector P B enforces the constraint in Eq. (26 1, 



namely, the occupation number nfg. of an a-spinon will be always equal to the occupation number of the 6-spinon in 
|RVB) with the opposite spin index a = —a at site I, if and only if there is a holon sitting at the same site I in |$^) 
(schematically it is represented by a dashed arrow in Fig. [2]). Here the phase factor e l ~ ls is introduced to simplify 
the formulation, whose physical meaning will become clear later. 

Corresponding to this exact mapping in the Hilbert space, the t-J Hamiltonian in the phase string representation 
can be further transformed as follows. 



2. The hopping term 



Rewrite the hopping term H t in Eq. ( |54| ) as 



(65) 



Further note that a shift A\- — > Af - — SAfj has to be made, when the 6-spinon Hilbert space is changed to a neutral 
spinon basis satisfying ^ ct n b icr = 1 instead of Eq. (51 1. Here 8A S ^ = | J2i^ij l^iQ) ~ ^j(0] (So- an \a) nl i IS obtained 
based on the definition (57). Then under the mapping (64) the hopping term is correspondingly changed to 

H t ^ PH t 

with 



) + h.c. 



(66) 



with P = PbP s where P s enforces the constraint in Eq. (15). 

In obtaining Eq. (66), the extra phases arising in Eq. (65) has been absorbed by making the following choice: 



(67) 



with utilizing Eq. (26). It is easy to verify that the phase factor e % ~ ls will serve as a 2D Jordan- Wigner operator to 
precisely make the a-spinons defined in Eq. (64) be have as fermions. Note that in the earlier approach |17j. without 
introducing e 1 " 13 to absorb the extra phase in Eq. (65), the backflow spinons are treated in a boson representation 
which is not as compact as in Eq. (66). 

Based on the precise mapping (64), one has 

|*> ^ Pb [|*.) ® |RVB) ® |$ h )] 

= P[l$„)®l* 6 )®l*h)] 

= P|*), (68) 
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where |$ a ) denotes the pure a-spinon state and |RVB) = P s \$i,). Generally \iff) here should be understood as a state 
expanded in terms of the direct product bases of the a-spinon, 6-spinon, and /i-holon. Hence in the hopping term 



(69) 



H t directly acts on the fractionalized state | "J) defined in Eq. ( 68 ) 



Corresponding to the mapping in the Hilbert space as given in Eq. ( 68 1 , the electron annihilation operator defined 
in Eq. ( 48 ) in the phase string representation may be reexpressed by 



Ci „m i— > p 



(70) 



which is obtained with using Eq. km, resulting in the phase shift field Q; = ($f - $"W2 defined in Eq. (KL2k. 

Similarly the spin oper ator s Si — in the phase string representation [cf. Eqs. (52 1 and (53)1 can be rewritten 



under the constraint Eq. (|51J) as S, = S|(l - nf) = S£ - n^Sf . Then under the mapping of Eq. ( |68| , it can be shown 
that 



tfS? I*) h> —PS? 



with using Eqs. (64 1 and (67), such that 



SiW^pfsJ + s?]!*), 

where the a-spinon spin operators S? are defined by 



and 



S 



a+ _ 



(71) 



(72) 



(73) 



(74) 



(we assume the anticommutating relation, e.g., a^a^ = — a^oL, between the a-spinons of opposite spins without loss 
of generality). 



3. The superexchange term 



For the superexchange term in Eq. (55), by introducing the factor (1 — nf)(l — nA to explicitly enforce the no 



double occupancy constraint (51) in the enlarged Hilbert space y&), one finds the following mapping 



H.m — ► -^p[(l-nf)(l-^)(A| J .) t A| 



where 



PHj\q> 



(75) 



(76) 



i?j in Eq. ( 76 ) can be further reexpressed as 



(*i> <u> 
+J^(S?-S!; + Sf -S!?) + JN h 



(77) 
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in which in obtaining the second term on the right-hand-side (rhs) of the first line, the relation 



A Q 



according to Eq. (64 1 is used, and in obtaining the second line the equality 



A s - = 



-Sj • S) + 1/4 



(78) 



(79) 



as well as Eq. ( 71 ) are utilized 



Therefore, by introducing a new kind of "backflow" a-spinon, the original fo-spinon in the phase string formalism will 
now always describe a "half-filled" neutral spin state constrained by Eq. (15). Consequently two distinct processes, 
involving spins in the hopping and superexchange terms of Eqs 
separately in terms of two kinds of spinons, as in Eqs. (66) and 
Simons gauge fields, Af, 



and (55), can be mathematically depicted 
In this new formalism, the mutual Chern- 



and A^j, are still defined by Eqs. (57) and (1581), but Afj is now always acting on a half- filling 
background. The phase difference in A\^ is then absorbed into the backflow spinon, and naturally turns the latter 
into a fermion. In the following, we will see that this kind of fractionalization description will be very important for 
properly constructing a saddle-point state at low doping. 



C. Mean- field scheme 



The t-J model has been first reformulated in the phase string representation in order to accurately keep track of its 
peculiar and unique sign structure and then expressed in a specific fractionalization formalism in terms of a neutral 
bosonic spinon, a backflow fermionic spinon, and a bosonic holon, in order to properly distinguish the microscopic 
hopping and superexchange processes in the restricted Hilbert space. Now one is ready to construct an effective 
theory/ground state based on the above new formulation. 

The Schrodinger equation H t -,j\^) — E\^) can be rewritten as 



P (H t -j - E) I*) = 



(80) 



by using H t -j\^f) i — > PH t _j\^} and |*) H> P|*), where 



H, 



t-J 



Ht+Hj 



(81) 



as defined in Eqs. (66) and (77). 

Based on Eq. (80), one may further make the following ansatz that 



H t -j\*)=E\V) 



(82) 



holds for both the ground state and low- lying e xcita tion states. Generally speaking, it is a sufficient but not a 
necessary condition for P\^f) determined by Eq. (82) to be an eigenstate of PH t -j- Later in Sec. Ill D we shall 



show that the Bogoliubov quasiparticle excitation is indeed an exception, which satisfies Eq. ( 80 ) but not Eq. ( 82 ) , 



emerging as a collective mode beyond the latter. But in the following we first focus on the solution of Eq. (82) and 
develop a generalized mean-field scheme. 

Such a mean-field theory will be underpinned by a gauge-invariant bosonic RVB order parameter 



A- 



(83) 



which was first introduced in Refs. jTOl HH] ■ Then Hj in Eq. ( 77 ) can be linearized in terms of the "mean-field" order 



parameter Ay, giving rise to an effective Hamiltonian H s in Eq. (34) with the order parameter (Af^) taken as 
s- wave- like: 



<:/ 

fA 



A s and J s = JA s /2. Self-consistently one always finds the mean-field 



b} a h, 



= 



NN 



(84) 



such that (A^) is the unique order parameter for the bosonic RVB state [101 HE]- Equation (84) implies that the 
bosonic RVB pairing is indeed incompatible with the hopping, in contrast to the fermionic RVB case, which has been 
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the basis for introducing the backflow spinon to facilitate the hopping process in the first place [T7]- It is noted that 
H s remains invariant, when (A^.j is changed from s-wave to d-wave-like, via a simple gauge transformation: 6j CT — > 
(-l)**&i ff so long as Eq. ((84]) holds. 



The kinetic energy of the a-spinons will arise from the hopping term H t in Eq. ( 66 1 , where a natural gauge- invariant 



decoupling gives ri se t o a pure holon term in Eq. (33) and the a-spinon hopping term in H a [i.e., the first term 
on the rhs of Eq. (35)]. The rest of terms in H a , including the pairing term for the a-spinons, all come from Hj in 



Eq. (77). The pairing term in Eq. (35) is obtained by a conventional mean- field decoupling: 



t 



h.c 



(85) 



term is incorporated into the first term in Eq. (35)], with a d-wave order parameter 



[a linearized e 

Afj = (<&„| Afj\& a ) — rjij |A"-|. Therefore, a "pseudogap" state characterized by the bosonic RVB pairing of Eq. 
which persists from half-filling to a finite doping, will be described by an effective Hamiltonian H e g = H s + Hh 
as given in Eq. d32h [Sec. II B2]. 



(1831, 



It is important to point out that H e g obtained based on Eq 
tonian. For instance, at the mean-field approximation, one would have th 



(|82|) is not strictly a conventional mean-field Hamil- 

**«\ = tK a and 



tE 



,t 



tn = t 



h\hje 



iA",, \ — 



= tH , but under the projection P, one will have 



t a ~ t at low doping. In principle, 
these parameters together with J s and J a , appearing in Eqs. (33), (34), and (35), are not determined by a standard 



self-consistent mean-field procedure. Instead, they will be c onsidered as the variational parameters, to be determined 
by minimizing the ground state energy of P^q) in Eq. ( 20 ) with regard to the exact H t -j. Such a variational scheme 
will decide the magnitudes and doping-dependences of these parameters, but their detailed values are not crucial to 
the qualitative physical consequences discussed in the present work. 



D. Quasiparticle as a collective excitation 



The above mean-field treatment of H-t-j leads to an effective Hamiltonian £f e g in Eq. (32 
fractionalized superconducting ground state (31 ). Based on H c g, two types of unconventiona 



, which determines the 
elementary excitations 

have been identified, i.e., the spin-rotons and the fermionic a-spinons, which have been discussed in Sec. II CI and 
C2, respectively. 

A quasiparticle mode, whi ch carries both charge and spin-1/2 quantum numbers, can be created by the electron 
operator defined in Eq. (70). In the superconducting phase, with the holon condensation \h\) ^ and the phase 



coherence ( e lUi ) ^ 0, an a-spinon excitation may directly appear in the single-particle spectral function according to 



Eq. ( 70 1 . But such a decomposition structure is only a part of the low-energy feature around the antinodal region as 
discussed in Sec. II C2 and C3. 

In the following, we demonstrate that a more conventional Bogoliubov quasiparticle will also emerge as an inde- 
pendent low-lying excitation in the superconducting state, in addition to the above non-BCS-type excitations. Being 
coherent around the nodal region, such a quasiparticle can be regarded as a bound state forming from the elementary 
building blocks of a holon, an a-spinon, and a nonlocal phase shift based on Eq. (70). Such a "collective" mode goes 



beyond the description of H c s in Eq. (32), and will be correctly described based on the original full t-J Hamiltonian. 



The hopping term in the latter will provide the necessary binding force for a stable nodal quasiparticle, which has 
been previously shown [2Pj 127] in the phase string representation by using the equation-of-motion method. Below we 
give a similar proof based on the present decomposition given in Eq. ( 70 ) . 
A quasiparticle excitation state may be generally constructed by 



|k^) qp = 



t 



(86) 



where the creation operator at is a linear combination of c_k-<r and c k . The quasiparticle spectrum is given by 



£ k = (ka\H t ^j\ka) qp -E G 



(87) 



using ako-l^c) — (assuming \^g) being normalized). 
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For the t-J model, one generally has[ 



j=NN(i) 



j=NN(i) 



and 



[Hj,c u 



J 



E (i 

j=NN(i) 



J I 



j=NN(i) 



By acting them on the ground state in Eq. (20) and using the d-wave order parameters: 
K = X) CT ( c iCT C jci') t° linearize the rhs of the equations, one finds 



(89) 



(ci-fCji) and 



[H t -j,c ia ] |* G > - \tc& C J' CT + ^ Cia l* G ^ ~ J E ^W ac ]-<j\^G) + decay term + scattering term (90) 



j=NN(i) J j=NN(i) 

where t c g = t(l + 5)/2 +JK/4 and /i is the chemical potential. 



If one can neglect the higher order terms on the rhs of Eq. ( 90 1 , including a "decay term" and a "scattering term" 



to be given below, then Eq. (|90j) reduces to a linear equation in the c-operators, which can be diagonalized via the 
following Bogoliubov transformation: 



where u£ + u£ = 1, with = 1 + (et — /x)/i?k, = 1 — (ek — M)/-Ek, and 2ukfk = Ak/-Ek- It leads to 



with a BCS like spectrum 



Here £k = — 2i e ff (cos fc^ + cos k v ) and 



£ k = V(£k-M) 2 + (A k 



A k = 2 J ^(cosa^ + cosg J/ )Ak+ q 



(91) 



(92) 



(93) 



(94) 



which can be easily shown to be d-wave if A^ c is a d-wave superconducting order parameter [57|. [Note that Ak 

is scaled by an additional factor \ho\ 2 oc S in Refs. [21] [27], which is actually an artifact in linearizing [2 7j the 
equation-of- motion.] 

By noting that the chemical potential /i is given by requiring (j2ka c ka c ^^ = 0- ~ <^)^j one finds the quasiparticle 

state c^I^g) behaves similar to a Bogoliubov quasiparticle in a d-wave BCS state, where it is built on a normal 
state with a large Fermi surface satisfying the Luttinger theorem as decided by £k = As a matter of fact, in the 
above equation-of-motion calculation, if, say, the next-nearest-neigbor hopping t' is added to -f/t-j, a corresponding 
change in the band structure of ek will take place to give rise to a modified Fermi surface just like in an ordinary 
band theory. By contrast, t 1 can be effectively renomalized to zero in Hh and H a due to the phase string effect, as to 
be shown elsewhere, such that it will not directly affect the spin-roton and a-spinon excitation spectra. This further 
illustrates the distinction between a fermionic a-spinon and a true Bogoliubov quasiparticle. 

Hence a Bogoliubov quasiparticle is always a coherent excitation at a sufficiently low energy as protected by the 
ODLRO ([6J). As pointed out above, this Bogoliubov quasiparticle can be viewed as a bound state in the fractionalized 
ground state, whose wave packet can eventually break down at a higher energy in the superconducting state where 
the phase coherence ^ is still maintained. To properly understand such a stability, one needs to further inspect the 
higher order terms in Eq. ( 90 1 . 

The "decay term" in Eq. ( 90 ) has the following leading contribution 



decay term 



(95) 
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after a linearization with using the order parameters: Hq, K a , A?-, and h^, as well as ^e lS1 y. Here the coefficients are 
given by O — 3tK a — JH , Pij — \tH a + ^JK a a, and Qij = | JA°-. Such a decay term represents the fractionalization 



tendency of a quasiparticle into an a-spinon, which is already indicated by the decomposition form ( 70 1 . It corresponds 
to a coherent term so long as the phase coherence (e 1 ^ ^ is maintained, in which the a-spinon appears at the antinodal 
region, and thus in the momentum space is distinguished from the Bogoliubov quasiparticle in the nodal region. Such 
a coherent term disappears above T c , even though the a-spinon can still remain coherent at T > T c . 
Furthermore the scattering term reads 



scattering term = 

j—NN (i) 



|*g> (96) 



which involves the scattering between the quasiparticle and the spin-roton excitations composed of the 6-spinons. 
By noting that a spin-roton excitation has a finite "resonance energy" E g in the superconducting state (cf. Sec. 



II CI), the scattering in Eq. (96) may be safely neglected if one only focuses on the nodal quasiparticles below 
this characteristic energy scale. On the other hand, at an energy scale higher than E g , a strong scattering between 
a Bogoliubov quasiparticle and the background neutral spin excitations is expected to dominate the single-particle 
spectral function in an appropriate momentum region, which will be important to understand the ARPES data, but 
is beyond the scope of the present paper. 

IV. CONCLUSION AND PERSPECTIVE 

The present mechanism for superconductivity resembles, by nature, what has been proposed by Anderson p] [28] that 
the spin RVB pairing is turned into the Cooper pairing upon doping. But the basic structure of the superconducting 
ground state, presented in Eq. is distinct from the original proposaljT] of the Gutzwiller-projected BCS state 
in Eq. ([!) by that the neutral RVB and Cooper channels remain clearly differentiated throughout the underdoped 
superconducting regime, as illustrated by Figs. [2] and [I] respectively. Such a ground state has intrinsically embedded 
three essential types of correlations: i.e., the AFM correlations in |RVB); the Cooper pairing of the doped holes; and 
a mutual influence/competition between the two channels via a mutual statistical phase in A^. 

Consequently the superconductivity arises as a self-organization from these correlations, rather than by default 
as in a Gutzwiller-projected BCS state at finite doping. Without holes, the long-range AFM order will always win 
in | RVB), which provides a highly accurate description of the ground state of the Heisenberg Hamiltonian. But a 
sufficient concentration of doped holes will eventually turn the antiferromagnetism in |RVB) into a true spin liquid 
state with short-range AFM correlations. And by doing so the Cooper paired holes can gain phase coherence to realize 
a high-T c superconductivity. 

The key underlying such peculiar structure in the ground state can be attributed to the singular phase shift 
introduced by the doped holes, which is incorporated in Eq. ^ by Ah- Such a "phase string effect" induced by the hole 
hopping is irreparable, representing the exact sign structure of the t-J model in the no double occupancy constrained 
Hilbert space, which appears as mutual statistical signs drastically different from the Fermi-Dirac statistical signs. 
Incorporating this new emergent statistics is thus essential in order to correctly understand the non Fermi liquid 
nature of doped Mott insulators. 

Correspondingly there are three distinctive types of elementary excitations. Firstly, the short-range RVB correla- 
tions in | RVB) are reflected by a finite energy gap E g opened up in a neutral spin excitation, called a spin-roton, 
which only reduces to the gapless spin wave in an AFM state with a long-range RVB pairing, e.g., at half- filling. Such 
a characteristic E g will protect the superconducting phase coherence from phase disordering effect, caused by A/j that 



closely monitors the neutral spin correlations in |RVB). In particular, it determines[25 a unique T c formula (42). 
Secondly, the fermionic excitations related to breaking up the Cooper pairs in Eq. ([3]) are called backflow spinons 
(a-spinons) as they carry well-defined spins and are therefore also contribute to spin excitations. Such backflow spinon 
excitations will constitute a lower branch of spin excitations, below E g , as the size of the Cooper pairing is usually 
larger than that of the RVB pairing. Thirdly, the a-spinon is distinct from a conventional Bogoliubov particle by a 

singular phase shift factor e li \ Although the phase coherence (e l ^j ^ in the superconducting phase can make the 

a-spinon show up in the single-particle channel around the antinodal regime, the local singular fluctuations of e lf2 will 
in addition induce a new collective mode as a bound state of the fractionalized particles around the nodal regime. It 
is nothing but a conventional Bogoliubov nodal quasiparticle. 

Some additional unconventional properties can be inferred in such a d-wave superconductor. In the spin channel, 
two branches of spin excitations near (ir, 7r), separated by a resonancelike energy scale E g , are responsible by the 
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spin rotons (upper branch) and a-spinons (lower branch), respectively. In particular, strong magnetic fields should 
affect the lower-branch first via the Zeeman effect, which can eventually destroy the pairing of the a-spinons and 
thus the superconducting order parameter, leading to small Fermi pockets formed by the fermionic a-spinons; On the 
other hand, in the single-particle channel, a dichotomy of Bogoliubov quasiparticle and the a-spinon coherent peaks 
can appear simultaneously in the superconducting phase in the nodal and antinodal regions, respectively. But the 
latter should disappear in the single-electron spectral function once the phase of the superconducting order parameter 
is thermally disordered by spin excitations in |RVB), where each excited spinon in |RVB) automatically induces a 
current vortex via e , forming a spinon- vortex and proliferating in a quantum vortex liquid state|24j: Furthermore, 
in the presence of weak magnetic fields, the Cooper pairing of doped holes provides a small phase stiffness p s (0) oc 6 
at T = 0. But p s will be thermally reduced by the Bogoliubov nodal quasiparticles which couple to the external 
electromagnetic fields with a full electric charge, leading to p s {T) — p s {0) — aT with a ~ 0(1). Obviously all of 
these will have strong experimental implications. But we shall further make the detailed comparison with the cuprate 
superconductors elsewhere. 

Finally, we remark that the present low-energy effective theory described by H c g in Eq. ( 32 ) resembles the so-called 
two-fluid models, which have been phenomenologically proposed for the cuprate superconductors [29, 30 and iron- 
based superconductors [31] based on different theoretical considerations. The main similarity with these approaches 
lies in that there exist both an itinerant BCS-type component (i.e., in the present case) and a localized spin liquid 
component (i.e., |RVB) here). The main distinction lies in how the local spin liquid component is mathematically 
characterized: the full description of it in the present approach is actually given by |RVB) ® \$h)> which is described by 
a mutual Chern-Simons gauge model [Eq. (40)] that respect spin rotation and time-reversal symmetries [19:. With the 
doping effect self-consistently incorporated, such a local spin component can naturally evolve from an AFM ordering 
state to a spin liquid state, and exhibit a multilevel pseudogap behavior in the underdoped regime [TO]. 
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Appendix A: A hidden ODLRO in ground state (Tsb 



The ground state ^ reduces to Eq. ([8| under the holon condensation condition ^ . If the total number of holons 
is not fixed, the ground state ([8| may be further expressed in a compact form 



\*a) =e^|RVB) 



Since the operator D will introduce a pair of holes, satisfying 



N h , D = 2D 



(Al) 



(A2) 



where N h = r$, by using N h (dY |RVB) = 2n |RVB), one finds N h \V G ) = 2D\^ G ) such that 



(E G |g|tt G ) = (N h ) SN 
<*g|*g> 2 " 2 



(A3) 



at a finite doping concentration S. The relation (A3) clearly indicates that the ground state (|8j) possesses a new 
ODLRO as given in Eq. ( 11 ) for a finite spatial separation of ij. 
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